1// noise in hopping conduction: Role of multi-site aggregates 
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We propose a mechanism for l//-type noise in hopping insulators based on the multi-electron 
charge redistribution within the specific aggregates of the localized states located in the vicinity of 
the critical resistors. We predict that the noise with 1 / /-type spectrum extends down to practically 
arbitrarily low frequencies. 
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A nature of the underlying mechanism of 1// low- 
frequency noise spectrum in hopping conductors is a 
part of the more general problem of the universal ori- 
gin of 1//, and had been an unresolved issue for 
decades i2iiiiSiSiL2i2iiS A general intrinsic reason for the 
1// spectrum is the exponentially wide distribution of 
the relevant relaxation times in the system. Thus the 
issue of the mechanism of 1// spectrum reduces to the 
question of the microscopic origin of such a distribution. 
The low temperature transport in doped semiconduc- 
tors occurs via electron hopping between the donors (or 
hole hopping between acceptors when considering p-type 
semiconductors). The hopping distance is determined by 
the balance between the tunneling probability over this 
distance and the probability of the thermal activation 
necessary to accomodate the difference between the ini- 
tial and final levels. Optimization of these probabilities 
cuts the effectively conducting band around the Fermi 
level. This hopping mechanism is called variable range 
hopping (VRH) transport. With the decreasing temper- 
ature this effective impurity band gets narrower and the 
electron has to tunnel over a larger distance to find the 
appropriate destination state, since the levels at the close 
by donors are separated by a large gap, see Ref. and 
references therein. In this regime the conducting paths 
form a dilute percolation cluster, and the system resis- 
tance is controlled by the very rare critical hopping re- 
sistors Re- The conductivity of the critical resistors are 
most susceptible to the dynamic fluctuations of charge 
distribution in a body of the conductor and they become 
the source of the current noise. 

The crucial question associated with the origin of the 
noise is whether the 1// law goes indefinitely to low fre- 
quencies and if not, what is the nature of its lower boimd. 
This problem was first addressed by Shklovskii& and Ko- 
gan&Shklovskiii who attributed the origin of the noise 
to the charge traps affecting the electron concentration 
on the percolation cluster that supports the current in 
the conductor. This model predicts the saturation of 
the noise spectrum at low frequencies. The reason for 
the saturation is that no trap can be perfectly isolated 
from the rest of the hopping sites since their spatial dis- 
tribution is Poissonian. Consequently, to secure a very 
large relaxation time one has to find a very rare trap 
specially isol ated from the other hopping sites. These 
ideas of Refs. were elaborated in the recent paper by 



ShklovskiiiS devoted to the temperature dependence of 
the noise spectrum. The alternative idea was put forward 
by Kogan^ who assumed that a hopping system is in fact 
the Coulomb glass and can therefore have arbitrarily low 
intrinsic frequencies due to transitions between the dif- 
ferent metastable "pseudo-ground" states separated by 
unlimitedly large barriers. Accordingly, there is no sat- 
uration in the noise spectrum down to extremely low 
frequencies. Unfortunately, no analytical description for 
this mechanism was proposed. 

Here we revisit the problem of the low-frequency noise 
in hopping conductors and develop the approach for 
calculating the noise spectrum based on the idea of 
multi-site bistable aggregates. While switching between 
their pseudo-ground states, the aggregates produce low- 
frequency electrical noise, which is "read-out" by the sites 
belonging to the backbone percolation cluster. The read- 
out mechanism suggested earlier— attributes the noise to 
fluctuations of the effective resistors forming bonds of the 
percolation cluster with correlation length C. The fluc- 
tuations of each resistor are due to the noise produced by 
the nearest bistable aggregate. Since the switching times 
of these aggregates are spread in an exponentially broad 
interval the noise spectrum turns out to be of 1// type. 

In the original approach^ the aggregates consisted of 
pairs of sites. However, the probability of finding such 
pairs - fiuctuators - is limited by the fact that such fluc- 
tuator be located in a "pore" free of other sites able to 
facilitate electron hopping^ Consequently, the maximal 
switching time are limited that leads to a cut-off in the 
noise spectrum Here we generalize the approach— includ- 
ing into the consideration multi-site aggregates that can 
assume two distinct charge configurations with the nearly 
same energy. In this respect, the aggregates that can be 
viewed as two-level systems. It is essential that transi- 
tions between these states occur via the multi-electron 
hops, either coherent or incoherent. The probabilities 
of such transitions arc strongly suppressed comparing to 
the case of isolated pairs, and this is the source of very 
large intrinsic switching times. Consequently, the cut-off 
frequency can be extremely low. The suggested model of 
multi-site bistable aggregates is a "bridge" between the 
ideas of small fluctuatorgSiiS and of pseudo-ground state 
of the whole system— 
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I. MODEL 

To quantify the concept of an aggregate let us start 
from the simplest one - a fluctuator consisting of the two 
trapping sites occupied by a single electron. Such pairs 
have two levels with small energy separation E <T and 
their relaxation time, Tp, is related to the tunneling be- 



tween the states; t 



e 2rp/a -^i^gj-g jy is gome (gen- 



erally temperature-dependent) pre-exponential factor, Vp 
is the spatial separation between the sites while a is the 
localization length. The values of t should be much larger 
than the typical hopping times, Th = I'Q^e^^'^/"- , along the 
current-carrying percolation cluster, therefore rp should 
be larger than the typical hopping distance in the clus- 
ter, r/i = a^. Here ^ is the connectivity parameter of 
the percolation theory entering the exponent of critical 
hopping resistor, Rc = Roe^, see Ref. fill and references 
therein. We consider an isolated pair, therefore there is 
no any other site in a close neighborhood (i. e. within the 
distance f <C r-p) from either of the centers, that could 
facilitate the processes with the rates 3> Tp"^. One can 
express this requirement as the inequality 



2r„ A 2r 
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where A is the relevant energy band width. To re- 
late this bandwidth to the distance r we estimate f as 
f = (3/47™)^/"^ where n is the concentration of the 
relevant centers. Since for a Coulomb-gap controlled 
system in the VRH regime n = (2/3)(KA/e^)^ we ob- 
tain A ~ (9/87r)^/^(e^/Kf). Optimizing the r.h.s. of 

Eq. 1^ with respect to f we find Tp < ThC^'-^^^K Here 
e = (To/T)i/2, To = e^na^ 

Now we consider the aggregates consisting of 4 hop- 
ping sites occupied by 2 electrons. These close to square 
arrays are the simplest next hierarchy level aggregates 
realizing the two level system (the three-site aggregates 
are reduced to pairs, since even a slight deviation from 
the perfect arrangement lifts the triple degeneracy and 
gives rise to an exponentially wide gap separating the 
split levels of the closest pair from the third level). A de- 
generate configurations requiring the least energy arise 
when two electrons are trapped at the diagonally oppo- 
site sites. Labelling the sites as 1, 2, 3 and 4, we say that 
the pairs are formed by the occupied sites 1 <-!■ 2, and 
3 4, respectively. For the configuration where sites 2 
and 4 are occupied one has 
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the configuration where sites 1 and 3 are occupied obeys: 
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Here et are the single-particle energies of the sites, r^j are 
the sites spacings, and /i is the chemical potential. The 
energy splitting between the two configurations is 

E = Ei-E2 = 2—[ +^2+£4-£i-£3. 4 

K \7'24 ri3 J 

The relaxation processes associated with the transitions 
between the two configurations with the almost equal en- 
ergies are the process that involve the simultaneous ex- 
change of electrons between the sites 1 — > 2 and 3 — * 4. 
All other processes increase the Coulomb energy and 
therefore are unlikely to occur. The switch can be ei- 
ther due to a multi-electron (a two-electron in our case) 
tunneling first suggested by Pollal*i^ or due to incoherent 
processes. 

As follows from the above mentioned correlation, the 
aggregate formed by the sites with large single-particle 
energies would have short intersite distances and, ac- 
cordingly, high rates of the multi-electron hops. Since 
the slow relaxing aggregates should have relatively large 
larger intersite separations, their single-particle energies 
should be relatively small. 

To proceed further it is convenient to map the aggre- 
gates onto a spin system. A single pair is mapped on 
the 1/2-spin (the spin direction is taken from the oc- 
cupied to the unoccupied site), the four-site aggregates 
arc represented as a pair of the interacting antiparallel 
spins. The interaction between the spins forming an ag- 
gregate is antiferromagnetic since the Coulomb repulsion 
pushes electrons belonging to the adjacent pairs apart. 
Adding more pairs we map larger aggregates to clusters 
of 3, 4, ...N spins. This model is a minimal model for the 
general Coulomb glass since we allow only for the transi- 
tions that are represented by flips of spins. In the absence 
of randomness the system is antiferromagnet, disorder 
may, in principal, drive it to the spin glass. 

The 2Af-site aggregates (or antiferromagnet A^-spin 
clusters) must obey restrictions similar to those given by 
Eqs. (O and l|2Jl. To estimate the density of states for the 
aggregates we shall take into account that in course of the 
reversal of the "spins" the energy of each site crosses the 
Fermi level. Thus we should compare the change of the 
site energy resulting from the aggregate rearrangement 
and the "static" scatter of the site energies. 

To do so, let us specify some energy ej within the 
Coulomb gap and consider the sites with energy band 
£i < £/. Let us consider the aggregates formed from 
these sites in the same way as it was considered above. 
For the site energy we have, cf. with Ref. 0, 



</5^ + V^, V, 



Here ipi is the potential induced on the site i by the back- 
ground charges not included in the aggregate; the sum 
in Vi is calculated over the sites forming the aggregate. 
Thus ifi is due to the disorder while Vi depends on the 
state of the aggregate. 
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To estimate Vi let us split it into the parts having 
different symmetry with respect to rearrangement of 
the aggregate between its two metastable states: Vi = 
(1/2)(1/+ + V~Si). Here s,; = ±1. The symmetric part 
can be treated as a contribution to the potential tpi . As- 
suming the potentials ipi and V^^ to be random and un- 
correlated, we write 



1 



(6) 



Since we are interested in the states within the Coulomb 
gap, the density of states for the single-particle energies 
Ei is given by the Efros-Shklovskii law, g oc ef. Thus the 
typical distance between the sites as 



g{e)de 



'1/3 



= (3/2)i/3eV 



(7) 



One expects that the potential formed on the site i by 
the other sites forming the aggregate can estimated as 

/ ^\ 1/2 

V± = [{^i ~ C±Sf where C± ^ 1 Biic numerical 

parameters depending on the aggregate geometry. Mak- 
ing use of Eq. © one has 

^ = e2(l-CV4-e/4). (8) 

Since the typical variation of the site energy, 5V = V- = 
C-Sf we obtain the following estimates 



C-£/, V+ = C+£f, ^ = Ef^l-e/4-Cl/4. (9) 



Note that the relations between V+ and (f> do not 
depend on provided we deal with the states within 
the Coulomb gap. One can expect that these ratios are 
of the order of unity. 

Basing on the discussion given above let us estimate 
the aggregates density of states. First, let us note that 
for all the sites belonging to the aggregate the inequality 
Vf > + 2(p,; must be met. Indeed, only such states 
cross the Fermi level The probability to form an aggre- 
gate depends on two parameters - the ratios C-/C+ a-^d 
C,^/2tf. Assuming that C+/C- 1 a-nd > one con- 
cludes that to form an aggregate one should deal with the 
sites with the potentials ipi less than their average scatter. 
Since the Coulomb gap exists only for the single-electron 
energies, Ei rather than for disorder-induced potential iy9, 
the distribution of the random potential ip , Vi^ip), is 
smooth for small (p and one can put V{(p) w l/2(^. Con- 
sequently, he relative weight of the aggregate sites is in 
this case given by 



1 

2^ 
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= A. 



These considerations are valid until A < 1. In this case 
relative weight of the aggregate formed by 2N sites is 



X2N ^ g-2Arini/A^ -jj^^g at A < 1 is the probability 
to form a metastable aggregate is exponentially small in 
terms of the number of the sites involved in the aggregate. 
This will be mainly discussed in what follows. 

However, there is no fundamental principle requiring 
A < lii^ This is why later we will also discuss the sce- 
narios where the multistable aggregates can be realized 
with a probability close to unity. 



II. THE DISTRIBUTION OF THE 
RELAXATION RATES. 

To calculate noise in the system one needs density of 
states, P{N,r, E), of finding the aggregate with N sites 
and distance r between the sites per unit energy interval, 
E. This density can be expressed in the form 



PiN,r,E) 



A 



N 
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(10) 



Here Nr"^ is the volume of the aggregate while \fNe^ j kt 
is the energy bandwidth for small energies. 

Let us estimate the distribution function. V(y)^ of the 
relaxation rates, ;/, for the aggregate rearrangement. We 
define the rate as 



V = vi^ max \ e 



(11) 



Here p = r/a, ^ = {Tes/TY^^. The first term in paren- 
theses describes formation of a "domain wall" in the ag- 
gregate and the second term corresponds to coherent tun- 
neling transitions leading to re-charging of all aggregate 
sites. 

To calculate the distribution of switching rates let us 
introduce the hmction L{N,p) as 

LiN,p) = -lnmax{e-^'''«'/2.7p^g-2Arp| ^^3) 

Thus L{N,p) ~ ln(i'o/i') = A^- Denoting pjv by the 
equality N^/^£_^ /2.7pN = 2NpN, 

PN = 0.43e/A^^/' , (13) 
one can express L{N, p) as 

L{N,p) = -lnmax{e-"-«6^'''«(''"/''),e-°-86^'''«''/''")} 
Let us now invert the above dependence. Since 

L{N,p)'. 



0.86iV5/'^e(p/pAr) , p:$>pN 

p:$> PN ■ 

the inverted hmction, p{L), has two branches: 
pjL) _ 

PN 





dp 


PN 


n, 
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dL 


Ln 1 
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(14) 



(15) 
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Here Ln = QMN^^^^, thus pn/Ln = 1/2N. Note that 
L{N,p) < Ljq. Now we are ready to calculate the distri- 
bution function 



r{u,E) = 47r y dNr'^drP{N,r,E)5[v-voe 



L{rM) 



dN p^dp P{N, p, E)5[p- p{L)] 



dp 



Substituting P{N, r, E) from Eq. (fTn|l wc obtain 



V{v) = 



1 27r-2.7 

V TesV^ 
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7V5/2 



M 



dL 



(16) 



where 7 = Inl/A, while N^iM) = {M/OM^^^^. Sub- 
stituting these values to Eq. (|16|l we get, 



1 1 4V^-2.7 
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(17) 



In this approximation the density of states V does not 
depend on the energy splitting E. The expression (I17f) is 

valid at -fN^ = [{'y^/y0.86^)\n{i^o/iy)]^^^ > 1, or 

H^oM>^h'^'. (18) 
The product jNc can be expressed as ln(z^o/i^)" where 
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From realistic estimates one can expect that the second 
factor does not significantly change within the experi- 
mentally accessible frequency interval. For example, as- 
suming i/Q = 10^° Hz we get M^^^ = 1.87 at = 1 Hz 
and Ai^^^ = 2.06 ai v — 10~^ Hz. Consequently, for rea- 
sonable experimental conditions the quantity Oi{v) can be 
replaced by a = Oi{v) where v is some characteristic fre- 
quency within the measurement interval. As follows from 
the above estimate, for any feasible frequency a ^ 1 and 
the noise spectrum is of the 1// type. 

As a result, the distribution of the relaxation rates can 
be expressed as 



C 



v^-^W{vo/v) Tes 



(19) 



where C ~ 12. 2/7 is a numerical factor. The temperature 
dependence a oc ^^^/^ cx T^/^ can be used for verification 
of the proposed mechanism. 



III. ESTIMATE OF THE NOISE INTENSITY 

Switching between the aggregate configurations leads 
to the change in the energies, Sef^, of the sites belong- 
ing to the percolation cluster and induces the current 



noise. If |fc-'^''| > T then the fluctuations of the hop- 
ping resistor, pi, are large, \5pi\ ~ pi and can change the 
percolation cluster structure (see Ref . 0) • 

To be specific we consider an A'^-spin aggregate coupled 
to the hopping resistor pi formed by the two hopping sites 
which we label by indices 1 and 2. The fluctuations of 
the resistance are estimated as 



-min(l,|fe,|/T) , 5e, 
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(20) 



If the distance between the resistor i and the nearest 
aggregate, rj, is much larger than the typical distance 
between the sites belonging to the hopping cluster, Vi ^ 
rh, the fluctuation in energy Ssi « {dSe^'^'^ / dri)rh where 
Se is the potential induced by the aggregate. The latter 
can be estimated assuming that the total dipole moment 
of the relevant aggregate is erc^/Nc where Tc = apj^^ — 
0.43aCA^c"^^*'. Consequently, = OASr^Nc^^^ and 



/a: 



-ES 



(21) 

As we have already mentioned the noise is formed at 
the exponentially rare critical resistors. Let us assume 
for a while that the fluctuations of their resistances are 
small, \5Ri\/ Ri <C 1. In a linear approximation the total 
resistance fluctuation can be written as 



5R 1 ^ 5Ri 



(22) 



Here M is the number if critical resistors in the sample 
which can be expressed through its volume, V, and the 
correlation length of the percolation cluster, C — = 
a^^, as A/" = \> I C? . Assuming the aggregates acting upon 
different critical resistor to be statistically independent 
we get. 



{6R,SR), 



{5R^,5R^ 



(23) 



Now let us estimate of 5Ri. One has in mind that the 
aggregates are also rarely distributed. This means that 
all the critical resistors have at most one neighboring 
aggregate and the fluctuation spectrum is: 



R? 



1 



Here Vi is the switching rate of the nearest aggregate 
with energy splitting Ei^ v^{ri) = min {5e{ri) /T)^ , 1 



the squared dimensionless coupling. Since 5e{ri) cx r~ 
there exists the specific distance, r^, at which the energy 
variation given by Eq. H21|l is equal to T: 



^^4/3^1/9^,, ^^1/3^1/9 



(24) 



At Ti ^ rx the interaction strength w^(ri) decays at least 
as cx r^^ ■ Consequently, only the nearest aggregate is 
important. 
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Replacing summation over the critical resistors (and 
their nearest aggregates) by averaging and integration 
over E we obtain 



Sic.) ^ 



Here 



2N 



Viuj). (25) 



Jj^l^^^ f^(r) d'^r, while the distribution 



function V{uj) is given by Eq. or Using Eq. lf?T|) 
we conclude that is of the order of unity, and 



Siio) 



AT 



(26) 



The factor ^/Af = ^C^/V = (aVV)C^ cx T-^^. Thus 
the the considered mechanism leads to a strongly de- 
creasing temperature dependence of the noise intensity. 
The suggested procedure to estimate the noise intensity 
is valid if each cell of the backbone cluster contains only 
one aggregate. This implies the limitation to the cluster 
size 



Ncri < C 



(27) 



This requirement hold for any realistic frequency since 
^ 3> 1. One can imagine another restriction relevant to 
the percolation mechanism behind the resistivity Accord- 
ing to the 1/uj spectrum, one expects that the magnitude 
of the fluctuations increases with an increase of the ob- 
servation time. If the relative resistance fluctuation for 
any relevant resistor, 



becomes comparable with unity, then one expects that 
the percolation cluster is completely reconstructed by the 
fluctuations. As a result, the fluctuations with lower 
frequencies will not be observed in the sample resis- 
tance. As follows from Eq. (|26(l . this requirement pro- 
vides the frequency limit for the validity of our calcu- 
lation, InivQ /lu) 3> VJ; which is automatically fulfilled 
for all the frequencies less than the typical hopping fre- 
quency, Vh = i^oe^^. 



IV. COULOMB GLASS SCENARIOS 

The above considerations demonstrate many-site mul- 
tistable configurations in the hopping insulator. The 
probability of finding such configurations depends on 
the interplay between the Coulomb interaction and the 
random potential produced by stray disorder, e. g. by 
charged acceptors. The question which arises here is 
whether there is a critical value of disorder discriminating 
between the cases of exponentially low, cx e^^^ , proba- 
bility of finding such configurations and the probability 
close to one. Earlier we relied upon rather strong disor- 
der assuming exponential decay of the probability with 
the size of the metastable aggregate. 



Now we will discuss the consequences of the intrin- 
sically correlated distribution of charges where there is 
no exponential decay of the probability with the expo- 
nent proportional to the number of sites in the aggre- 
gate. Namely, let us assume that there is a critical value 
7c of the parameter 7 = Inl/A such as that the sys- 
tems with 7 < 7c allow to form the arbitrarily large 
metastable configurations. The independent support of 
this idea (although again for model systems) is given by 
theoretical considerations given in Ref. Namely, it 
is demonstrated that the sites within a Coulomb gap 
demonstrate replica-symmetry breaking which can be re- 
lated to a presence of extensive number of metastable 
states within a thermodynamic limit 00. We also 

note that multi-stable character of the ground states was 
also demonstrated by several numerical simulations&i^ 

In principle, one can consider the following 3 scenarios: 

(i) Strong disorder, ^ > "fc- Only rare multistable ag- 
gregates are available due to the effect of disorder, 
the probability exponentially decreasing with an in- 
crease of the aggregate size. 

(ii) Weak disorder, 7 < 7c. The multistable configu- 
rations are inherent for the system, which in this 
case can be called the Coulomb glass. The results 
depends on the typical distance of charge transfer: 

a. the charge transfer within the metastable ag- 
gregates is restricted by neighboring sites. 
Correspondingly, the interactions are domi- 
nated by the short range dipole-dipole forces 
("dipolar glass"), and the short range ordering 
in the sites occupation numbers still exists. In 
this case one expects that the local glassy dy- 
namics is not significantly affected by coupling 
with the remote regions. Namely, the dynam- 
ics of a critical resistor is dominated by the 
nearest aggregate with a specific relaxation 
rate. One can expect that in this regime the 
system would have the 1 /to spectrum down to 
any practically achievable frequency, and the 
exponent a — 0. 

b. The charge transfer at large distances is im- 
portant. In this case the local short-order 
configuration of the sites occupation numbers 
can depend on the state of the remote regions 
( "large-scale Coulomb glass" ) , and result 
a hierarchical dynamics becomes possible. In 
this case the low-frequency noise cannot be 
regarded as a superposition of statistically- 
independent Poissonian telegraph-like fluctua- 
tions occurring in different parts of the sample 
and acting on different backbone resistors. 
Because of hierarchical dynamics of the aggre- 
gates one can expect that fluctuation of each 
site energy possesses 1 / /-type spectrum. To 
the best of our knowledge, this situation has 
not been analyzed yet in a convincing way. 
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The difference between the "dipolar" and "large- 
scale" Coulomb glasses can be in principle revealed 
experimentally by analyzing local occupation fluc- 
tuations by means of scanning tunnelling spec- 
troscopy. For the scenario (iia) the telegraph noise 
with a given relaxation time is expected, while for 
(iib) a noise with a complex spectrum should be ob- 
served provided the large enough observation times 
are possible. 

For both scenarios, the 1// noise spectrum should hold 
down to arbitrary small (from practical point of view) 
frequencies. 

Since the characteristic switching times for large ag- 
gregates can be too long to be observable, at realistic 
noise frequencies one deals only with switchings of rel- 
atively small metastable aggregates. To consider these 
aggregates as compact and statistically independent one 
has to ensure that the aggregate's switching between 
the metastable states does not cause repopulation of the 
states outside the aggregate. Consequently, the signs 
of their energies (counted from the chemical potential) 
should not change. In addition, the Coulomb interactions 
with the sites outside the aggregate should not affect the 
aggregate dynamics. Both contributions are proportional 
to the surface of the aggregate, i. e. cx N"^/^ . Repeat- 
ing the analysis of Sec. [H] and ascribing the probability 
factor ^ '' to allow for the surface contribution one 
concludes that S{uj) (X uj~^ and a = 0. 

The controversial experimental results on the tempera- 
ture dependence of the flicker noise in the Coulomb glass, 
see e. g. Refs. Il6l can, in principle, be of fundamental 



nature. Indeed, in the previous discussion we did not con- 
sider the time fluctuations of the single-particle energies 
of the sites forming the aggregates. These fluctuations 
are due to correlated electronic hops, which can be more 
important for the noise than for the stationary trans- 
port. An attempt to consider the role of the fluctuations 
of site energies in hopping transport was made in Ref . IT^ 
Somewhat later an important role of strong fluctuations 
in the Coulomb glass was also emphasized in Ref. We 
feel that correlated hopping may indeed introduce a fre- 
quency cut-off to the noise spectrum; this cut-off can be 
considered as a hallmark of correlated hopping. This may 
be the case in the materials where experimental values 
of To are significantly smaller than those conventionally 
expectediii 

Our conclusion regarding the absence of the frequency 
lower cut-off agrees favorably with available experimen- 
tal results. The temperature dependence is more prob- 
lematic since different temperature dependences were ob- 
served, see, e. g., analysis in Ref. fioi 
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